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ABSTRACT 

We carry out a comprehensive theoretical examination of the relationship between the 
spatial distribution of optical transients and the properties of their progenitor stars. By 
constructing analytic models of star-forming galaxies and the evolution of stellar populations 
within them, we are able to place constraints on candidate progenitors for core-collapse 
supernovae (SNe), long-duration gamma ray bursts, and supernovae la. In particular we 
first construct models of spiral galaxies that reproduce observations of core-collapse SNe, 
and we use these models to constrain the minimum mass for SNe Ic progenitors to ~ 25 Mq . 
Secondly, we lay out the parameters of a dwarf irregular galaxy model, which we use to show 
that the progenitors of long-duration gamma-ray bursts are likely to have masses above 
» 43 M0 . Finally, we introduce a new method for constraining the time scale associated 
with SNe la and apply it to our spiral galaxy models to show how observations can better 
be analyzed to discriminate between the leading progenitor models for these objects. 
Subject headings: methods: analytical - gamma rays: bursts ~ stars: evolution - supernovae: 
general 



1. INTRODUCTION 

Transients are any events that exhibit a brighten- 
ing, dimming, or otherwise noticeable change within 
a finite and usually short lifetime. They are one of 
the most important observable objects in astronomy 
as they represent a major change of state, which of- 
ten exerts an enormous impact on its surroundings 
{e.g. McKee & Ostriker 1977; Dekel & Silk 1986; 
Gehrz et aL 1998; Matteucci & Recchi 2001). Fur- 
thermore, these events inform a great deal on the 
physics and evolution of the progenitor stars that 
cause them. 

Yet, despite their enormous importance, the pro- 
genitors of many transient types remain poorly con- 
strained. One would like to be able to make predic- 
tions about these events and watch them unfold in 
real time, but transients, as their name suggests, are 
fleeting and fickle - many observed only accidentally 
and often only partially. In particular, in the case of 
supernovae (SNe) and gamma ray bursts (GRBs), 
observers must rely either on extremely demanding 
direct methods or approximate indirect methods to 
constrain their progenitor systems. 

In the direct approach, one seeks to identify the 
progenitor star responsible for a particular transient 
on a case by case basis. This, in turn, requires 
one of two things: (i) the occurrence of a rare, 
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nearby event, as in the case of SN1987A (Gilmozzi 
et al. 1987; Kirshner et al. 1987) in which the sur- 
prising properties of the progenitor {e.g. Arnett et 
al. 1989) coupled with neutrino detections (Bionta 
et al. 1987; Hirata et aZ. 1987; Mayle et al. 1987) 
led to important advances in our understanding of 
stellar physics; or (ii) a laborious search through 
archival images to extract the one in « 10^^ stars 
per year that results in a supernova {e.g. Man- 
nucci et a/. 2005). Even with this second, "brute- 
force" approach, a good deal of serendipity is still 
necessary. Thus, until recently, only one object 
(1993J; Aldering et aZ. 1994), had been identified 
from this approach; the situation improving to a 
handful of objects with the advent of painstaking 
searches though pre-explosion Hubble Space Tele- 
scope (HST) images (Barth et al. 1996; Van Dyk 
et al. 2002; Smartt et al. 2004; Maund et al. 2005; 
Hendry et al.2006; Li et al.2006; Gal- Yam et 
al. 2007). 

A second approach relies on indirect constraints 
and uses the properties of transient host galaxies, 
or the locations of transient within these hosts, to 
derive progenitor properties. Studies of this type 
include measurements of the rates of supernovae 
{e.g. Pain et al. 1996; Cappellaro et al. 1999; Tonry 
et aZ. 2003; Blanc et a/. 2004; Dahlen et aZ. 2004; 
Sullivan et al. 2006; Barris & Tonry 2006; NeiU et 
al. 2006; Botticella et al. 2008; Dahlen et al. 2008) 
and gamma ray bursts {e.g. Daigne et aL2006; 
Charay et al. 2007; Liang et al. 2007; Guetta & Delia 
Valle 2007; Kistler et al. 2008), constraints from the 
metal content of galaxy clusters {e.g. Renzini 1999; 
Loewenstein 2001; Scannapieco et aZ. 2003; Tozzi et 
al. 2003; Baumgartner et al. 2005) and high-redshift 
quasars (Barth et aZ. 2003; Dietrich et a/. 2003), 
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studies of the distribution of SNe relative to spi- 
ral arms (Maza & van den Bergh 1976; Delia Valle 
& Livio 1994; Bartunov et aL 1994; McMillan & 
CiarduUo 1996; Petrosian et a/. 2005), and studies 
that relate the presence and properties of transients 
to the environmental properties of their hosts {e.g. 
Oemler & Tinsley 1979; Wang et al. 1997; Hamuy 
et al. 1996; Fruchter et al. 1999; Howell 2001; Le 
Floc'h 2003; Christensen et al. 2004; van den Bergh 
et al. 2005; Ostlin et al. 2008). 

Recently, Fruchter et al. (2006, hereafter F06) de- 
veloped a new such indirect observational method, 
building on the framework first laid out by Baade 
(1944) and elaborated on for GRBs by Bloom et 
al. (2002), and they successfully applied it to dif- 
ferentiate between the progenitors of long-duration 
gamma-ray bursts and supernovae. The method 
involves observing the spatial locations of GRBs 
and type-II SNe in their host galaxies and calculat- 
ing the fraction of the total host galaxy light con- 
tained in pixels fainter than these locations. After 
several such observations, a pattern emerges, and 
F06 demonstrated that long-duration GRBs have a 
higher propensity to cluster in the brighter regions 
of a galaxy than type II SNe. Kelly et al. (2007) ex- 
panded upon this analysis and applied it to type la 
and Ic observations, to find that type-Ic SNe have 
similar environments to GRBs. 

Theoretical work has followed to use the F06 
method to put constraints on which stars might con- 
tribute to GRB observations (Larsson et aZ. 2007). 
In this paper, we show how the F06 method, com- 
bined with simple analytic models can be used to 
derive strong constraints on the progenitors of a 
wide range of transients. In particular, we simu- 
late model galaxies that act as hosts for three vari- 
eties of transients: core-collapse SNe (type II's and 
Ic's), type la SNe, and long-duration gamma-ray 
bursts. Using star-formation rates and stellar evo- 
lution profiles, we construct simulated galaxies with 
all the features of a dynamical system of stellar 
births and deaths and demonstrate that the tran- 
sients are primarily determined by their progenitor 
life-times, and thus progenitors can be constrained 
by the locations of the transients. By simulating the 
locations of stellar deaths for massive stars that re- 
sult in core-collapse SNe and long-duration GRBs, 
we are able to place constraints on the progenitor 
masses. And by simulating the location of low mass 
stars that result in white dwarfs and type la SNe, 
we are able to constrain the time scale for these 
transient events, as well as suggest ways that the 
F06 method can be modified to provide further con- 
straints. 

The structure of this paper is as follows. In §2, 
we lay out our approach for simulating transients in 
model spiral galaxies, carry out tests on the resul- 
tant data, and apply it to derive mass constraints 
for core-collapse SN of various types. In §3, we lay 
out the parameters and approach for a simulated 
dirr galaxy, establish and test the criteria it must 



follow, and apply the method to derive mass con- 
straints for long-duration gamma-ray bursts. In §4, 
we return to our spiral galaxy model, employ analyt- 
ical models for the distributions of type la SNe in 
order to constrain the associated dynamical time- 
scale, and suggest a modified method to improve 
these constraints. Our conclusions are summarized 
in §5. 

2. CORE-COLLAPSE SNE 

Type lb, Ic and type II supernovae are all core- 
collapse SNe. That is, they are the result of a star 
whose core has collapsed to form a neutron star or 
black hole [e.g. Woosley et al. 1995). In such ob- 
jects, the core no longer produces enough energy 
to counteract its gravitational collapse, and it re- 
leases tremendous amounts of neutrinos, eventually 
resulting in a cataclysmic outflow of the surround- 
ing envelope {e.g. Burrows et al. 1995; Mezzacappa 
2005). The prevailing model for lb and Ic SNe is 
that of an asymptotic giant branch star that has 
been stripped of most of its envelope by a stellar 
wind before undergoing core-collapse (Wheeler & 
Harkness 1990; Woosley et al. 1995; Yoon & Langer 
2005; Woosley & Heger 2006). On the other hand 
type II's retain their gaseous envelopes during the 
core-collapse and thus feature strong hydrogen ab- 
sorption lines flanked by P-Cyg absorption troughs. 

Kelly et al. (2007) compared observed type II and 
type Ic counts against the integrated luminosity 
of the host galaxy in the B-band. In order inter- 
pret these observations, we first construct a suitable 
model for the integrated light of a test galaxy and 
compute the surface brightness in a coordinate grid 
(r, 9) for some time t. With this model in hand, we 
calculate the SNe rate for each pixel in the (r, 9) 
grid and compare the result against the integrated 
light. 

Most core-collapse SNe are found in late-type, spi- 
ral galaxies of middling to high metallicity, and to 
this end, we construct a model of a spiral galaxy 
with solar metallicity, Fe/H = 0.02. Note that Kelly 
et al. (2007) constrained the effect of the bulge on 
measuring fractional intensity in the B-band to less 
than 2%, and so for our purposes, we need only 
model the disk component. 

2.1. Approach 

To construct such a model, we start with an ex- 
pression for the star formation rate that takes the 
form of a rotating, density wave. In our model, 
we approximate the density wave as a Dirac delta 
function with a strength given by some radial profile 
0(r), 

OO 

t{r,9,t) = (j){r)^5{9-mT-npt), (1) 

ri=0 

where E is the star-formation rate per unit area, 
which, like (/)(r), has units of mass per unit time 
per unit area. In ([1]), f2p is the rotational velocity 
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of the star-forming wave, and the sum over delta 
functions ensures a rotating wave with two arms tt 
radians apart. 

To calculate the integrated luminosity at some 
time after the formation of these stars, we em- 
ploy the population-synthesis code, GALAXEV 
(Bruzual & Chariot 2003), which evolves a sample 
population of stars and outputs time-dependent lu- 
minosities, L^{t), in a chosen band, v. From this, 
we obtain an expression for the surface brightness 
of a pixel (r, 9) at a time t, 

E^(r,6l,t)= / L^{t-t')±{r,0,t')dt' 



\np\-\{2) 



The image created by this simple expression repre- 
sents a snapshot at a time t of the host galaxy and 
can be compiled for any band. 

To compute the SNe rate using the same expres- 
sion for SFR, we return to GALAXEV, which also 
calculates the main-sequence cutoff mass, i.e. the 
stellar mass that leaves the main sequence over time 
via the Padova 2000 isochrone (Girardi et al. 2000). 
A simple inversion then yields the main-sequence 
turn-off time with mass, t„is{M)^ which provides a 
more accurate version of the approximate expres- 

(\ —2.5 
years {e.g. Hansen et 

al. 2004) . Since this time interval is approximately 
90% of the total hfetime of a star (Tinsley 1980), 
this relation places the deaths of large stars near 
the locations of their formation, while lower mass 
stars die further away from their birthplaces. 

When stars form in a molecular cloud, large stars 
are less likely to form than smaller stars, and so, for 
a given population, a histogram of stellar masses 
will be heavily weighted toward the low-mass end 
(Salpeter 1955; Chabrier 2003). Observationally, 
one can thus construct an initial mass function 
(IMF) which relates the numbers of stars with their 
mass for a population, and throughout this paper 
we use the estimate by Chabrier (2003), which is 
also used by GALAXEV: 

(log Af-log _A/e) 




dhiM 



M < Mq 
M > Mq 



(3) 

where a = 0.69, and = 0.08 M©. At around Mc, 
the gravitational collapse of the star-forming molec- 
ular cloud is not as able to counteract the internal 
pressure, preventing cloud cores that are smaller 
than the Jeans mass from forming into stars and 
flattening the IMF. 
For high mass stars, the supernova rate, ^'(i) is, 

-m = ^ = imuoM)r'^, (4) 

where Mms(i) is the stellar mass corresponding 
to a main-sequence turn-off time t. We can at- 
tenuate this expression with Heaviside functions 



e(M - M- ) and e(M - A/+) to set the lower and 
upper mass limits for a given SNe type. The final 
expression is then 



^'(t, Af-,Af+) l.l[M{Q.9t)] 



-2.3 



dM 

~dr 



X [e (M - A/+ ) - e {M - M-)] (5) 

We convolve this expression with the SFR in the 
same fashion as with the integrated light to produce 
a per-pixel supernova rate at a time t for the model 
galaxy. 
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I],„(r,0,i,M-,Af+) = 
d0 [ '^{t~t')t{r,9,t')dt' (6) 



^ 1 1.1, 



x[e(Af-A/+)-e(Af-Af-)] (7) 

Returning to the expression for the radial profile, it 
is simple enough to begin with an exponential in r. 



(t){r) = So exp 



— r 



(8) 



where tq is the scale radius and Eg sets the star 
formation amplitude. 

The pattern speed, 17p, of a spiral galaxy is de- 
pendent on the epicyclic frequency, k(t'), and in the 
frame of the stars, a representative model is 
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1 

71 



(9) 



vo/r r > Ikpc 
wo/lkpc r < Ikpc. 



This velocity will have a singularity at 



ro 1 



1 

71 



3.414ri 



0, 



(10) 



where the pattern does not move with respect to the 
stars. This is called the corotation resonance orbit 
(Goldreich & Tremaine 1981; Binney & Tremaine 
1987) and any simulation would have to terminate 
at radii below this singularity, which fortunately is 
far enough away from ro to not affect the outcome. 

Finally, the surface brightnesses per pixel in the 
images obtained from ^ and (O are sorted into 
two columns according to the brightest to lowest 
pixel in the integrated light image. The output of 
each integral is a 2D linear grid in (r, 9) , the effect 
being that the pixels are keystone-shaped, and so 
each pixel must be weighted by r when calculating 
the fraction of light dimmer than each pixel as is 
done in the F06-style plot. 

2.2. Tests 

To understand the behavior of this plot, it is nec- 
essary to isolate the two mechanisms that place 
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SNe co-spatially in pixels of a certain brightness: 
the natural dimming due to the radial profile and 
the location of star formation as opposed to stellar 
death. Within a dynamical system of star formation 
that depends only on an angular coordinate such as 
a grand-design spiral galaxy, constraining observa- 
tions to an annulus removes the ambiguity intro- 
duced by the radial profile. This ambiguity can be 
described as introducing a second measurement of 
SNe intensity at a particular brightness that is not 
due to stellar evolution, but simply to densities of 
stellar populations. In an annulus, which fixes the 
model to a single radius, irrespective of the radial 
profile, the coordinate is analogous to a time co- 
ordinate, and features evident along 9 are due only 
to stellar evolution. Unraveling the annulus, then, 
is akin to looking through geologic strata, and the 
locations of the SNe measurements can be seen as 
the age at which they occurred. 

If the model is further constrained to count only 
those SNe that are due to stars above a threshold 
mass, a pattern begins to emerge whereby the higher 
the minimum mass cutoff, the more concentrated 
the SNe profile in the brightest regions of the galaxy. 

One of the features evident in Figure [T] is the cut- 
off brightness for each of the higher mass thresholds 
where SNe corresponding to stars larger than the 
minimum mass are no longer observed. This tracks 
well with stellar age, and the brightness associated 
with each cutoff are a good indication of the speci- 
fied minimum progenitor mass. 

Another feature that must be accounted for is the 
nonlinear climb in the SNe profile with brightness. 
This is actually due to a second ambiguity in the 
brightnesses of pixels in the galaxy; for many pix- 
els of a certain brightness for each radius, there are 
two corresponding ages which can be on the order 
of 10-100 Myr apart. The effect of this is to intro- 
duce gaps in the comparison of SNe intensities to 
the integrated luminosity of the host galaxy. When 
plotted in a F06 fashion, these gaps result in non- 
linearities in the SN curve. The effect is small but 
non-zero and must also contribute to the integrated 
F06 profile. Figure [2] demonstrates this ambiguity 
in the B-band. 

Restoring the radial profile essentially smears the 
behavior shown in Figure [T] across the dimmer re- 
gions of the galaxy as it reintroduces the ambiguity 
of SNe measurements at dimmer locations coming 
from stellar evolution and, critically, from regions of 
the galaxy at large radius that simply contain fewer 
stars. 

Finally, a behavior that is most evident when con- 
straining observations to an annulus is the propen- 
sity for the mass-brightness relations to break down 
at longer wavelengths. At near- infrared wave- 
lengths, brightness is less strongly associated with 
star formation, and therefore the relationship be- 
tween a star's mass and its location within a galaxy 
is not as clearly evident. In these idealized models, 
we neglect the effects of dust, which will complicate 



the results to the extent that obscuration is biased 
to regions of ongoing star formation. The behavior 
of a sample of stars over time for several bands is 
plotted in Figure [3l 

2.3. Results 

Our first simulation, using solar mctalicities, 
counts SNe from all stars above 8Afg, which cor- 
responds roughly to the minimum mass for a core- 
collapse, type II SN. The total number of SNe per 
pixel can be thought of as a supernova intensity, 
with every SN counting equally toward the total in- 
tensity, i.e. the intrinsic luminosity of a particular 
supernova is not a factor. 

When comparing the B-band, integrated light 
profile of the simulated galaxy to a SN intensity plot 
for the same parameters of galaxy size and metallic- 
ity, as shown in Figures |4] and [5l it is immediately 
evident that the final locations of all stars larger 
than 8Mq produce a similar profile to that of the 
integrated B-band light of the entire galaxy. The 
majority of the galaxy's B-band intensity is emit- 
ted by the large stars that have relatively short 
lifetimes, approximately 35Myr for 8Mq, as com- 
pared to a typical galactic rotation time on the or- 
der of 300Myr. This places their deaths near to 
their birthplaces and hence to most of the B-band 
emission that occurred throughout their lifetimes. 

In both Figs. |4] and [5l the intensity peaks at or 
near the crest of the star-forming wave, where large 
stars are created, burn brightly in the B-band, and 
die rapidly. As we adjust our simulation to the 
redder end of the spectrum, this prominent feature 
which gives spiral galaxies their arms becomes less 
and less apparent, until in the K-band, the galaxy 
looks more like a featureless, exponential disk than 
a grand-design spiral. This also follows observa- 
tions that spiral arms become less pronounced as 
one moves towards longer wavelengths (Zaritsky et 
al. 1993; Rix & Rieke 1993). 

Returning to the B-band simulation for M > 
SMq , an F06-style plot of the pixel intensities (Fig- 
ure [8]) shows a correlation between type II SNe in- 
tensity and B-band light which was demonstrated 
observationally in FOG. If we adjust in (O to 
estimate the type Ic minimum progenitor mass as 
25Mq, the arm features in the SN intensity plot 
become much more constricted as compared to the 
B-band light. 

A plot of these two profiles yields a very similar 
behavior to that observed in FOG, namely that the 
ordinary type II core-collapse SNe show a nearly 1:1 
relationship with the B-band light. However, our 
simulation also demonstrates that observed type Ic 
SNe show a strong correlation to the spatial distri- 
butions of stars with masses above « 25 Mq. This 
supports the idea that type Ic SNe are the transient 
events associated with the deaths of large stars. 

A K-S test of the resultant data, comparing dif- 
fering mass cutoffs to the observations described in 
Kelly et al. (2007) and shown in Figure [9l confirms 
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Fig. 1. — A F06-style plot for B-band annuli with variable cutoff masses. The thick solid line shows a one-to-one correspondence 
between light and SNe. The thin solid lines correspond to SNe with various limits in M~, ranging from 4 Mq (an unphysical 
example shown for illustration purposes) up to 50 Mq, all with upper limits of 100 Mq. 



Time [Myr] 

100 200 



300 



0.8 - 



J 0.6 - 




0.2 [ 



e 

Fig. 2. — A B-band luminosity profile for an annulus. Luminosity is normalized by its maximum. 



that 25 M© is the most probable candidate for the 
minimum mass responsible for type Ic SNe. Sullivan 
et al. (2006) suggest that 15% of all core-collapse 
SNe take the form of type Ic's, and so, if 



/g °" M-^dM 



^ 0.15, 



(11) 



then Af < 30 M© is the largest minimum pro- 



genitor mass for type Ic's that could yield a 15% 
SNe rate, given the Kroupa IMF (Kroupa 2001). 
The gray area in Figure [S] indicates disallowed min- 
imum masses for type Ic's based on this analysis. 

In the type II case, on the other hand, the best fit 
is obtained when one imposes both a lower mass cut- 
off and an upper mass cutoff of « 20 M© as shown 
in the type II curves in Figure This fact is borne 
out in Figure [TUl which plots K-S probabilities for 




various ranges of masses and demonstrates a greater 
than 85% probability of a maximum mass of 20 Mq 
for type II SNe. Figure [10] also demonstrates the 
lack of any compelling upper mass cutoff for type Ic 
SNe, and so a similar maximum mass is not imposed 
for type Ic's. 

2.4. Observational Requirements 

Two important observational requirements exist 
for the analysis laid out here. The first is a min- 



imum image resolution that provides enough data 
for performing an accurate analysis, and the second 
is a maximum inclination that allows one to accu- 
rately distinguish between different regions within 
the galaxies. 

In the extreme case of too little resolution, one 
might imagine a host image constrained to a single 
pixel. In this case the transient location would not 
exhibit any dependence on stellar mass, as it would 
always occur at the same position. To test the ef- 



Fig. 5. — SNe intensity for M > 8Mq , corresponding to tlie B-band liglit distribution of tiie model spiral galaxy. 




feet of image resolution on our results, we compare 
results for the mass cutoffs found in §2.3 at vari- 
ous image resolutions, from 100 pixels, up to 10^. 
In this analysis, absolute pixel resolution and ef- 
fective resolution are not distinguished so as not to 
obscure the physical insight with observational com- 
plications. Figure [TT] demonstrates that the overall 
nature of the plot does not change once the pixel 
density goes above a minimum threshold of 500 pix- 
els. When we repeat the K-S test for the mass cut- 



off found in §2.3 using only 500 pixels as opposed to 
lO'^, the confidence level changes by only 2%. Be- 
low the critical threshold of about 10^-10^ pixels, 
however, there simply are not enough data to illus- 
trate a convincing relationship for the locations of 
transients. 

Inclination plays a complicated role within the 
F06 analysis. At the extreme end, with an edge-on 
spiral galaxy, the locations of transients can only be 
said to be near or far from the centers of their hosts, 



Fig. 7.— SNe intensity for M > 25M 
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Fig. 8. — F06-style plot for core-collapse SNe in spiral model. Observed type lb, Ic and II SNe from Kelley et al. 2007 are 
plotted for comparison. The diagonal line follows a 1:1 relationship between fractions of SNe and fractions of total light. 



rather than correlated with star-formation. At an 
intermediate inchnation angle, i, for a given reso- 
lution, the vertical extent of the galaxy contracts 
by cos(z), decreasing the total area, and so for each 
location within the host galaxy, the corresponding 
fraction of light increases. The global increase in 
the surface brightness does not play a role in the 
final FOG analysis, but the decreasing area has the 
effect of concentrating transients in smaller regions. 
Figure fT2l demonstrates this behavior for the type 



II and Ic SNe cases. For i < 45°, the effect of in- 
clination is marginal, but for larger values of i, the 
FOG profiles rapidly diverge from that of a face-on 
galaxy. Kelly (2007) specifically excludes "edge-on" 
galaxies, placing the upper limit on their sample 
host inclinations i < 70°. Excluding any sampling 
bias, this places the average inclination at i « 35°, 
which is below the threshold where such an inclina- 
tion could be distinguished from face-on in the FOG 
analysis. Moreover, because high inclinations con- 
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Fig. 9. — K-S test for type Ic and type II SN minimum progenitor mass. A better fit to the data is achieved is we impose 
a 20 Mq upper cutoff on the type II simulations. In the Ic case, no upper mass cutoff is imposed, but the gray area indicates 
disallowed minimum masses above which there would be too few Ic SNe to be consistent with the observed rates. 
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Fig. 10. — 2D K-S plots for type II (left panel) and Ic SNe (right panel) with contours at the 85% confidence level. While 
the left plots favor a finite upper mass limit for Type-II SNe, no upper mass cutoff is required to match observations of type Ic 
SNe. The shaded region below 30 in the right panel corresponds to disallowed minimum type Ic progenitor masses. 



centrate transients in brighter regions for reasons 
not related to star formation or stellar evolution, 
galaxies with high inclinations {i > 45°) are not el- 
igible for the F06 analysis. For example, an FOG 
analysis comparing observed type II SNe to a simu- 
lated host with i = 60° would indicate a minimum 
progenitor mass of « ISM©, which is not consistent 
with accepted models of stellar physics. 

A further complication that inclination could 
bring is that of the propensity for dust lanes to ob- 
scure the brightest regions of spiral galaxies at high 
inclination. Because we ignore dust and have pre- 
sented a case for avoiding high inclination spirals, 
this effect is not treated here. 

3. GAMMA RAY BURSTS 



Gamma-ray bursts are the brightest observable 
objects in the universe, far outshining their host 
galaxies, and they come in two flavors: short- 
duration bursts and long-duration bursts (eg., Kou- 
veliotou et aL 1993). The favored model for long- 
duration GRBs, the so-called collapsar model, is 
that of a very massive star whose core collapses 
to form a black hole (MacFadyen et al. 1999; Mac- 
Fadyen et al. 2001). As the stellar envelope is con- 
sumed by the newly- formed black hole, tremendous 
energy is released in the form of gamma rays along 
the axis of rotation, and it is this beaming that 
gives GRBs their brightness when pointed directly 
at Earth (Rhoads 1997). The light-curves and spec- 
tra of the afterglows of long-duration GRBs occa- 
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sionally match those of Ic's in their latter stages 
(Rhoads 1997; Galama et al. 1998; Delia Valle et 
aZ. 2003; Hjorth et a/. 2003; Malesani et 2004; 
Campana et al. 2006; Plan et al. 2006; Kelly et 
al. 2007). 

Short-duration GRBs, on the other hand, are less 

luminous than their long-duration counterparts and 
may or may not be the result of a coUimated beam- 
ing effect. The model for producing these events 
has not yet been settled, though possible candidates 
include neutron star mergers (Eichler et al. 1989; 
Narayan et al. 1992: Ruffert & Janka 1999) and 
magnctar crust failure (Hurley et al. 2005). For the 
purposes of this paper, only long-duration gamma- 
ray bursts will be considered and henceforth, simply 
called GRBs for brevity. 



3.1. Approach 

To accurately model the light distribution of GRB 
host galaxies, we must construct a reasonable ap- 
proximation of a low-metallicity, dwarf irregular 

galaxy that exhibits pseudo-random clustering of 
star-forming regions. The necessary properties of 
our model galaxy are: 

• The average radial brightness must follow an 
exponential profile (e.g., Barazza et aL2006). 

• Starbursts should occur at random inter- 
vals for random initial masses, but the time- 
averaged star formation should yield an ex- 
ponential profile in total stellar mass (e.g., 
Barazza et a/. 2006). 
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• For simplicity, each starburst in a particular 
location should be independent of starburst 
history; that is, the propensity for a starburst 
to occur in any location should not decrease 
if that location has already experienced active 
star formation. 

• Starbursts should follow the observed initial 
cluster mass function in their composition 
(Oey et al. 2003; Dowell et al. 2008). 

• Old stellar populations should be diffuse 
and less concentrated than young populations 
(Parodi et al. 2003). 

To build a model that fits these criteria, we first 
fix the total star formation as a function of radius 
with an exponential profile 



S(r) = /3 exp 



— r 
ro 



(12) 



where /3 is a tunable parameter that sets the total 
star formation amplitude and rg is a scale parameter 
which sets the half-light radius by the relation Rm ~ 
1.68ro. Next we fit the initial cluster mass function, 
which describes how many OB stars, Nqb, per unit 
mass are contained within a cluster, as 



dN 



dlnNoB 



■(xN, 



OB ■ 



(13) 



This can be reproduced by choosing each cluster to 
contain 

NoB{e) ^ {I ~ e)i^) , (14) 

where e is a random variable ranging from to 1 
and a = 1.88 is taken to match the distribution of 
OB associations in dwarf irregulars (Oey et al. 2003; 
Doweh et al. 2008). 

For the last criterion, that old stellar popula- 
tions will become unbound after « 30Myr as the 
gas is ionized and leaves the cluster (Blitz & Shu 
1990; Matzner 2002; Krumholz et aL2006), result- 
ing in stellar diffusion away from the cluster cen- 
ter, we first calculate the dynamical time for one 
of these clusters to completely diffuse. Given that 
the stars in a cluster will have some initial, pe- 
culiar velocity, Av, we can estimate the diffusion 
time in the tangential direction as ~ {nr/Av), 
where r is the radial coordinate of the star clus- 
ter within the galaxy. For our model, we use a 
value for Av on the order of lOkm/s (Hunter et 
aL2005). The extent to which the cluster has dif- 
fused if it is younger than td.e is simply given by 
A9 = {Av/Vc)[itci - SOMyr) /r], where t^i is the age 
of the cluster and Vc is the circular velocity of the 
stars at the galaxy's edge. 

For radial diffusion, the maximum extent to which 
the cluster can diffuse can be derived from a simple 
energy argument given that a dwarf irregular pos- 
sesses a solid body potential (Binney & Tremaine 
1987) with stellar velocities v{r) = rcu, where oj = 
Vc/R, with R being the radius of the galaxy and 



Vc = v{R). If the /'-coordinate were decomposed 
into X and y coordinates, dy/dt — rujcos{ujt), and 
with a maximal kick in the positive or negative y- 
direction, (dy/dt)' — {rui ± Av)cos{Lut) and y — 
(r ± Av / uj) sm{Lut) . Therefore, Ar^ax — (Av/lo), 
regardless of the radial coordinate of the cluster. 
The diffusion time in the radial direction, then, is 
td,r = Ar/Av = 1/uj and is the shorter of the 
two diffusion times, and for our purposes, can be 
taken to be instantaneous after the cluster becomes 
unbound. During diffusion, the stars will be ar- 
ranged according to a Maxwell-Boltzmann distribu- 
tion with 



dn 



Ar 



r {RAv/VcY 



■ exp 



Ar 



RAv/Vc 



dn 
dA0 



AO 



(tsAv/r)- 



■ exp 



A0 



tsAv/r 



(15) 



, (16) 



where tg here refers to the time since the cluster has 
become unbound. 

For each radius, we compile a list of starburst 
events at random times with a random 9 coordi- 
nate until the total number of stars reaches that re- 
quired by (fT2)) for that radius over the lifetime of the 
galaxy. This list is then combined with GALAXEV 
and equations (fT5|) and (fT6|) when rendering the 
galaxy to incorporate the evolution of the star pop- 
ulations formed in each starburst, the time- variable 
brightnesses for each OB association, and the dy- 
namical diffusion of each cluster. The result is a 
model that fits our established criteria for a simpli- 
fied dirr galaxy, four examples of which are shown 
in Figure [T^ 

3.2. Tests 

Before we apply our model to derive F06-style 
plots, we must first ensure that our simulated galax- 
ies match star-formation profiles for observed dIrr 
GRB hosts, which are generally brighter and have 
more active star formation than ordinary dirrs. We 
accomplish this by adjusting /3 in to place the 
resultant log(SFR/kpc^) vs. Mb within the GRB 
host dataset, shown in Figure [U with GEMS and 
SDSS observations plotted for comparison (Barazza 
et al. 2006; F06). The SFR for the simulated galax- 
ies is a calculation from the resultant galaxy's to- 
tal stellar mass output over its lifetime rather than 
the input SFR. The SFR with error bars for the 
observed GRB hosts are range estimates given the 
radii and brightnesses of the host galaxies provided 
by F06, and assuming constant star formation for 
between 10 Myr and 1 Gyr. 

An important test our simulated galaxies must 
also pass is that of the relation between the scale 
lengths of the integrated light and of the star for- 
mation rate. If the star formation rate remains con- 
stant in time, as it does in our model, small stars 
formed at the edges of the galaxy, where star for- 
mation is rare, will persist in the integrated light 
spectrum for all bands, and this effectively stretches 
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Fig. 13. — Sample dirr galaxy simulations of B-band luminosity with diameters of 6kpc. 
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Fig. 14. — Properties of simulated high- and low-redshift dlrrs, 
dlrrs from GEMS and SDSS (Barazza et a/. 2006). 

the scale length for the integrated light beyond that 
of the star formation rate. The effect will obviously 
be most pronounced in the long-wavelength band- 
passes, but will also affect the short- wavelength 
bands as a portion of the output spectrum of old, 
red stars is also in these bluer bands. Parodi 
et al. (2003) constrains the relation between the 
SFR and B-band scale lengths for nearby dlrrs to 
'ro,sfr ~ (0.86 ± 0.06)ro,B. As this is a local con- 
straint it must be compared to dlrrs with somewhat 



as compared to the GRB host dataset from F06, and nearby 

lower SFRs, and thus we need to lower the SFR of 
galaxies as indicated in Figure [TH Measurements 
of simulated galaxies with star formation rates and 
brightnesses that correspond to nearby dlrrs yields 
''o.s/r = (0.90 ± 0.07)ro._B which is within the error 
of the Parodi et al. (2003) results. 

3.3. Results 

As with the type Ic SNe in the spiral model, 
long-duration GRBs show a strong correlation to 
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Fig. 15. — F06-style plot for long-duration GRBs in dwarf irregular model. Observed GRBs and type II SNe are plotted for 
comparison (F06; Kelly et al. 2007). 



high-mass stars when added to an F06-style plot. 
Figure [15] also repeats the earlier observation that 
type II SNe progenitors are strongly correlated with 
the B-band light, even in dwarf irregulars. A K- 
S test, shown in Figure [161 reveals that a mini- 
mum mass of « 43 Mq is the most probable can- 
didate for GRB progenitors. This value is in agree- 
ment with the findings of Larsson et al. (2007) and 
supports the conclusion that GRB progenitor stars 
must be massive enough to form black holes as op- 
posed to white dwarfs or neutron stars (MacFadyen 
et aL1999; MacFadyen et aZ. 2001). One result of 
particular importance is that though the F06 tracks 
for GRBs and type Ic's show a characteristic similar- 
ity, the masses associated with their progenitors are 
quite different. This suggests a dependence on the 
morphological properties of the host galaxy for the 
F06 method. For the type II SNe probability dis- 
tribution, again, an upper mass cutoff of 20 Mq is 
imposed, yielding a greater probability for « 7.7 M© 
being the minimum mass necessary for a type II SN. 

Note that the underlying physics in this com- 
parison is very similar to that in our spiral-wave 
modeling. In the case of the spiral model, the ro- 
tating density wave dynamics was prime motivator 
for concentrating high-mass, B-band emitters in the 
brightest parts of the galaxy. In the dwarf irregu- 
lar case, the dynamics of stellar clusters play this 
role, concentrating young, blue stars nearest the lo- 
cations of their birth and hence near the brightest 
portions of their host galaxies. 

4. TYPE lA SNE 

Lastly, we apply our approach to type la SNe, the 
so-called standard- candle events by which distances 
to host galaxies can be measured (Colgate 1979; 



Branch & Tammann 1992; Phillips 1993). These are 
thought to be the result of a white-dwarf accreting 
mass above the Chandrasekhar limit and collaps- 
ing into a neutron star either via stable accretion 
from a binary (Nomoto 1982; Starrfield 1994; Web- 
bink 1994) or a collision with another white dwarf 
(Woosley & Weaver 1994). Like type Ic SNe, type 
la's are often found in late- type, spiral galaxies, and 
we focus purely on the populations found in such 
hosts. In modeling the locations of these transient 
events, we must take into account the different pro- 
cesses that lead to type la's as compared to core- 
collapse SNe. Specifically, the incidence of type la's 
is not directly related to progenitor mass. 

4.1. Approach 

Type la's have been observed to have two com- 
ponents: a larger component that is dependent on 
the star formation rate and the time since star- 
formation, and the smaller component that is SFR- 
independent(Mannucci et al. 2005; Scannapieco & 
Bildsten 2005): 



SNRujt) 
(lOOyr)-i 



A 



Mr. 



B 



A'hit-r) 



, (17) 



where is the total stellar mass and M*(t — r) 
is the star formation rate at some previous time, 
t — T. Scannapieco & Bildsten (2005) presented the 
formalism for constraining the values of A and B, 
and Sullivan et al. (2006) find 5.3 ± 1.1 x lO'^^ and 
3.9 ± 0.7 X lO"'', respectively. 

This relation establishes all the necessary infor- 
mation for populating la events in our model galaxy 
without making any assumptions about progenitor 
masses. Whereas in the core-collapse model, stellar 
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Fig. 16. — K-S test for GRB and type II SNe progenitor masses. A 20 Mq upper mass cutoff is imposed on the type II 
simulation as in §2, while the GRB simulation has no upper mass cutoff. The observed, most probable minimum mass for 
GRBs, f5i 43 Mq, agress with the findings of Larsson et al. (2007). 

evolution and lifetimes were the dominant factors 
in placing the transient events in particular spatial 
locations, in the la model, the locations of the B 
component transients are controlled by a character- 
istic delay-time, r, {e.g. Madau et al. 1998; Gal- Yam 
& Maoz 2004) in the perturbation after star forma- 
tion. 

When plotting the la distributions in an F06 fash- 
ion for the B-band, the effect of changing the delay 
time for the B component is not immediately no- 
ticeable, as shown in Figure[T7l To better constrain 
the delay time of the B component, we suggest a 
modification to the FOG method. In this case, the 
image of each host galaxy is normalized by its over- 
all radial, exponential profile, which emphasizes the 
spiral density wave distribution with respect to the 
overall stellar distribution. That is for each galaxy 
we take 

t,{r,e,t) = j:,{r,9,t)/cP{r), (18) 

where (f){r) is the average radial profile of the galaxy. 
Similarly, each of the type la transient event images 
is normalized by this same exponential profile, with 



EsN(r,0,t) = SsN(r,0,t)Mr). 



(19) 



Each of these normalized distributions is then 
used to calculate the normalized light in pixels 
fainter than each galaxy, that is with central pixels 
artificially dimmed. Finally, these distributions are 
combined into F06-style plot using the normalized 
number of la's, that is with la's occurring near the 
outskirts of disk galaxies, contributing more than 
those near the crowded centers of these galaxies. 
This simple procedure draws out the spiral-wave 
structure, resulting in large differences between dif- 
ferent delay-time distributions. 



Because this method relies on transient locations 
being dependent on the ordered dynamical proper- 
ties of the host in conjunction with the delay time, it 
is most useful for applying to spiral galaxies rather 
than irregulars. Moreover, since ellipticals do not 
experience any star formation, they are also not el- 
igible for this analysis, or else r would have to be 
inordinately large. 

4.2. Results 

In Figure [191 we show the results of carrying out 
this modified FOB analysis on our model disk galax- 
ies. Given that regardless of the candidate pro- 
genitor model for type la's (binary accretion, white 
dwarf mergers, etc.), a WD star must have first been 
produced, the minimum value for r is necessarily the 
shortest lifetime of a white dwarf progenitor star. 
Core-collapse occurs for stars of mass greater than 
8Mq which corresponds to a lifetime of approxi- 
mately 35 Myr at solar metallicity, and so r > 35 
Myr. This minimum value is about 10% of the dy- 
namical time of the model galaxy, and thus we may 
vary r by as much as an order of magnitude while 
still retaining the usefulness of the FOG analysis. 

What is evident from the modified FOG plot is 
that with an increasing delay in the onset of the 
B component, the la profile crosses the integrated 
light curve at increasingly dimmer fractions of the 
total light. Figure [T51 demonstrates why this is so. 
As the delay time for the dynamical component in- 
creases, the peak output of the SNe intensity shifts 
to angular coordinates correspond to ever decreas- 
ing luminosities in the integrated light, causing the 
crossover in the modified FOG plot at lower fractions 
of total light. 
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Fig. 17. — F06-style plot for la SNe in spiral model with several delay times. Observed type la SNe are plotted for comparison 
(Kelley et al. 2007). 




Fig. 18. — A B-band luminosity profile for a spiral galaxy and a SNe profile for t = 35 Myr, respectively. The bottom pair 
have been normalized by an exponential to remove the non-dynamical component of the la SNR. 

This behavior, when matched against observa- 
tions that employ our modified analysis, should 
serve well to constrain the time delay between star 
formation and la transient events, thus placing lim- 
its on what kinds of progenitor models are viable 
explanations for type la SNe. The accretion model 
requires only that one white dwarf has formed in bi- 
nary with an evolved star and would be consistent 
with a value of r « 35 Myr, while the white dwarf 



merger model introduces a degeneracy and would 
likely push the value of r to longer time-scales. Fur- 
thermore, this modified procedure can be simply ap- 
plied to existing data sets, an application which we 
are actively pursuing. 

5. CONCLUSIONS 

The FOG method is a valuable tool in the search 
for transient progenitors. By relying on simple ob- 
servations, the method remains agnostic to the de- 




Fig. 19. — Modified F06-styIe plot for la SNe in spiral model, in which each pixel is normalized by </>(r). 



tails of the transient events, such as spectroscopic 
signature or peak intensity, and instead places con- 
cise and testable constraints on their spatial distri- 
butions. 

In this paper, we have demonstrated first, that 
models of spiral galaxies illustrate the dynamics 
that lead to the observed spatial distributions of 
core-collapse SNe. In a spiral galaxy, the rotat- 
ing star-formation wave is the primary dynamical 
component for setting the spatial distributions of 
SNe with varying progenitor masses. The distances 
that stars are able to travel from this wave during 
their lifetime determines how F06 profiles for cer- 
tain transient events will present themselves. Since 
there is a direct relationship between stellar lifetime 
and mass, our analysis has constrained type Ic SNe 
only to stars above w 25 Mj) and has supported es- 
timates for the minimum mass of type II SNe of 
«8M0. 

In the dwarf irregular model, the role of the star- 
forming wave is played by the star-formation pro- 
file which, in conjunction with the cluster diffusion, 
places young stars in the brightest portions of the 
galaxy. By employing GALAXEV results in combi- 
nation with observed dependencies on SFR and the 
IMF, our dwarf irregular model fulfills our initial cri- 
teria for being a representative counterpart of real 
galaxies as well as adhering to the observed inte- 
grated light profiles of dwarf irregulars. In combin- 
ing the F06 method with the results of this model, 
we are able to place constraints on the progeni- 
tor masses of GRBs to stars with masses above 
« 43 Mq . Our dirr models also confirm a mini- 
mum progenitor mass for type II SNe of « 7.7 M©. 
In future work, we will improve our simulation by 
adding a reasonable amount of patchy absorption in 



the dwarf and spiral galaxies. With too much dust, 
the GRB results would not conform to observational 
constraints, and so, these models might also serve 
to constrain the effects of dust in GRB hosts. 

For type la SNe, the delay time associated with 
the dynamical component of the two-channel la in- 
tensity profile is the primary factor in determin- 
ing their spatial distribution. However, the effect 
of changing the delay time for the dynamical com- 
ponent is not evident in a standard F06 plot. For 
this reason, we have suggested a modification to this 
analysis that normalizes each pixel by overall radial 
profile of the host. Applying this method to our sim- 
ulated galaxies demonstrates a clear dependence of 
the resultant profile on the delay time of the dynam- 
ical component. Future application of this method 
to existing observations of type la SNe will serve 
to place concise constraints on the delay time, sup- 
porting or refuting candidate progenitor models for 
SNe la. 

For each of the objects we have considered in this 
paper, we have demonstrated that galactic dynam- 
ics play an important role in the distributions of 
stars at the moment of their deaths. Furthermore, 
while not considered here, this techniques is likely to 
be applicable to other types of objects such as clas- 
sical novae, Wolf-Rayct stars, and planetary nebu- 
lae. From the observed distributions of transients, 
we are able to constrain the masses and properties 
of the progenitors from which they are formed. As 
transients flicker in and out of existence, the physics 
of stars is written across the night sky. 

We are grateful to Lars Bildsten, Bruce Doak, 
Brian Gleim, and Sangeeta Malhotra for their many 
useful comments and suggestions. We are also 
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